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DISCUSSION OF HYDRAULIC FUNDAMENTALS OF CLOSED 
CONDUIT SPILLWAYS 
PROCEEDINGS-SEPARATE NO. 354 


GEARY M. ALLEN, Jr.,! J.M. ASCE.—The writer agrees that the 
phenomenon of flow in a closed conduit system presented by the author 
may appear novel to many engineers; however, virtually the same flow 
condition has been confronted by hudraulic engineers engaged in the 
design of flood control outlets through large dams for many years. 

Of necessity, conduits constructed through the upper levels of ma- 
sonry dams are laid on a slope which becomes tangent to the face of 
the overflow section at the conduit’s exit. In most cases this design 
requirement results in a conduit laid on a steep slope which could pro- 
duce a flow condition through the conduit similar to that which the 
author describes. 

Of particular importance to the designer is the location of the hy- 
draulic control of the steep conduit with respect to reservoir stage. 
Generally, the location of the control with respect to reservoir stage 
is estimated by assuming the control at the outlet end of the conduit 
(full flow condition), computing the rating curve of the conduit for this 
condition and comparing this curve with the part full rating curve. 

The pressure or hydraulic grade line throughout the length of the 
steep conduit and especially immediately below the operating gates is 
of significance because as the author states, if the friction grade line 
falls below the conduit (conduit centerline in the case of fig 6a) pres- 
sures less than atmospheric will occur with resultant cavitation tenden- 
cies. It must be remembered that the hydraulic grade line represents 
the average pressures and not the maximum or minimum pressures that 
may occur, consequently, care should be exercised in the selection of . 
allowable average pressure reductions. 

The author is to be complimented on his description of the actual 
phenomenon, and on his presentation of a more general method for the 
determination of the pressures within the conduit. 

WEN-HSIUNG LI,? A.M. ASCE.— The phenomenon of pipes flowing 
full at steep slopes regardless of the outlet condition has been observed 
independently by the writer and by other investigators® in connection 
with studies of flow through culvert-pipes with rounded entrance. It is 
of great interest to the writer that this phenomenon has also been ob- 
served by the author on pipes with sharp-cornered entrance. 


1. Hydr. Engr., Corps of Engrs., Albuquerque, N. Mex. 

2. Asst. Prof. of Civ. Eng., The Johns Hopkins University, Baltimore, 
Md. 

3. Straub, L. G., Anderson, A.G. and Bowers, C.E., “Importance of 
Inlet Design on Culvert Capacity,” Research Report No. 15-B, High- 
way Research Board, 1953. 
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The advantage of having the pipe flowing full at steep slopes can be 
brought out forcefully by a graph similar to Figure a. This graph has 
been plotted for the determination of culvert-pipe size for known values 
of allowable depth of ponding h above the invert at the entrance, the dis- 
charge Q, the slope of the invert s, and the Manning’s roughness coef- 


shl/3 D 
ficient n of the pipe. When “a is greater than about 100 and h is 


from approximately 0.65 to 0.87, the flow in the pipe is pulsating be- 
tween full flow and part-full flow. The values of shown in the 
P 


graph are the minimum values observed. For pipes flowing full at slope 
s other than the friction slope, the values of 7% depend also on the 
length of the pipe. The values shown in the graph serve to indicate only 
the general trend of variation. As can be seen from the graph, larger 
pipe does not necessarily give larger discharge with given values of h, 
s andn. Under certain conditions, the same capacity can be obtained 
with pipes greatly = in size. For example, if the design condi- 

sh 

n2 

or a pipe with D = 0.6h may be used. The larger pipe will flow partly 
full while the smaller pipe will flow full with subatmospheric pressure 
in the barrel. Pipe-culverts with sharp-cornered entrance have not 
been observed to flow full under practical conditions. For this case, 
the capacity always decreases with the pipe size, and in the example 
above, the larger pipe must be used. Obviously, it is of great practical 
interest to find the general criteria for maintaining full flow in a closed 
conduit at steep slope. 

In studies of the hydraulics of closed conduits, the hydraulic grade 
line is generally plotted at a distance above the center line of the con- 
duit equal to the pressure head at the center line. This method is cor- 
rect only when the pressure is hydrostatically distributed across the 
cross-section. At the free outlet of a pipe flowing full, the pressure is 
not hydrostatically distributed due to the slight curvature of the sur- 
faces of the jet. However, the pressure is hydrostatically distributed 
at a very short distance upstream from the outlet. At this point, the 
pressure head at the center line is equal to D/2. Thus the hydraulic 
grade line should meet the crown of the conduit at the outlet as shown 
in Figure 6. The available head across the structure is thus equal to 
the difference in elevation between the free-water surface above the 
entrance and the crown of the conduit at the free outlet. Thus the value 
of z in Eq. 1 should be equal to the difference in elevation between the 
crest of the weir and the crown of the conduit at the free outlet. This 
correction may be very small for the example given by the author. 
However, this correction of D/2 in the available head may be extremely 
important when D/2 is not small compared with the available head, as 
in the case of culverts for highway drainage. 

The author has mentioned model testing to determine whether cavita- 
tion would occur in a given structure. It is of practical interest to have 
a method of predicting the occurrence of cavitation without resort to 


tions are such that = 800 and rue = 3.0, a pipe with D = 1.3h 
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model tests. The minimum pressure occurs at the vena contracta near 
the entrance. Since the loss of head upstream of this section is negli- 
gible, the minimum pressure can be easily shown to be (see Figure b) 


2 
(1 + K, (s - s,)L 


If c is used to denote the coefficient of contraction at the vena contracta, 


a@ can be shown to be z . This equation awaits verification by 
observation. 

ARTHUR L. COLLINS,’ A.M. ASCE.—In the Separate 354, and with 
reference to the “CONCLUSIONS”, second line, the word “may” could 
imply that the conduit could flow full or partially full without submer- 
gence. The writer prefers to argue that the conduit may not safely 
maintain a high efficiency unless it is submerged, and offers the siphon 
outlet as an example. 

In the northern California delta, between Sacramento and Stockton, 
there are about 500,000 acres of reclaimed submerged tidal lands. 
Large levees surround man-made islands with areas of a few hundred 
to several thousand acres. First the levee is constructed from mate- 
rial borrowed from outside the area, thus leaving a deep channel. Inside 
the area the excess rain and seepage is brought to a suitable location by 
a large drainage canal, where it is pumped over the levee. 

A typical modern pump installation, of which there are perhaps 100 
with discharge pipes of 18 to 30 inches, consists of vertical pump in- 
stalled at the terminal of the drain canal at the levee. A total of 200 or 
300 feet of pipe extends from the pump to the base of the levee, thence 
up the levee slope, thence nearly level through the levee where it again 
descends until the end can be submerged in the waste canal. 

The crest of the siphon will be about 20 feet above the surface in the 
drainage canal, and 15 feet higher than the waste canal, thus leaving a 
net lift of about 5 feet for the pump. It is a simple matter to anticipate 
the hydraulic gradient in the pipe system and the head available from 
the pump. However, water gas always forms from the high rotation of 
the pump impeller, and accumulates in the crest of the siphon. Also, a 
reduction of the pressure at the siphon crest causes the gas to expand. 
It eventually stops the operation of the siphon unless the system is so 
designed that the gas capsule is broken into small volumes and flows 
down the pipe with the water. Expensive mechanical installations have 
been provided in some cases to take care of this trouble, when the 
proper tilt to the pipe and a deflector in the pipe would have solved the 
problem in an inexpensive and simple way. 

In the same area there are hundreds of smaller sizes of pipes, 4 
inches to 12 inches in diameter, which are true gravity siphons, and 
which siphon the water from the waste canal back over the levee for 
irrigation purposes. In observing these gravity siphons and pump si- 
phons over a series of years, one will find conditions such as pipe ends 
rusted off, or extra low water, which will expose the pipe ends. When 


1. Consulting Engineer, Berkeley, California. 
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this occurs there is but one chance in a hundred that the pump siphons 
will not break and cause an enormous overload on the pump motor, due 
to an increase in the lift of from 6 or 7 feet to 20 feet. There was one 

exception where a 20 inch siphon did operate without a break when not 

submerged. 

What happens is that immediately after the seal is broken at the pipe 
edge, a gob of air enters the pipe as a wedge which widens and travels 
up the pipe until the down pipe is partially filled. 

On the other hand, the small gravity siphons are often observed to be 
discharging full and without submergence. Here again air will collect 
from the water and accumulate at the crest of the siphon and eventually 
block its operation if the air is not carried along with the water. 

The simplest example of a full pipe discharge is seen in the use of 
a hose where gasoline is siphoned out of a barrel. Also it is common 
practice to use a 2 or 3 inch bent pipe to siphon water from an irrigation 
lateral to the land level without requiring submergence. 

As stated above, the break away or hydraulic drop at the pipe end 
permits the air to wedge its way up the pipe. Its occurrence depends 
upon the shape, size, slopes, velocities and other factors. 

Reference is now made to Fig. 1, (b). Using this as an example it 
appears that the outlet of the conduit may be the equivalent of the out- 
fall of a submerged pipe end. The water in the outlet basin is not free 
to fall away, and expose the end of the conduit, and the air wedge cannot 
become effective. Such would not be the condition should the pipe or 
conduit end as in (a) of Fig. 2. 

Can it be that the prototype on a small scale tells only half the 
story? In the pump siphon an unnecessary loss of head of only a few 
inches may amount of several hundred dollars in an increase of annual 
power cost. It is possible that if the gravity operated conduit was re- 
quired to operate at its maximum efficiency, the practice of submerg- 
ing the conduit end would become the rule. 
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DISCUSSION OF PRESSURE SURGE CONTROL AT TRACY 
PUMPING PLANT 
PROCEEDINGS-SEPARATE NO. 361 


C. E. WITHERS.-—In discussing the pump characteristics of the 
Tracy pumps the author states that the pump performance data neces- 
sary for investigating the transient hydraulic conditions in the zone of 
“energy dissipation” and the zone of “turbine operation” were not avail- 
able from the pump manufacturer and had to be estimated from the 
curves of another double volute pump. Due to the expense involved in 
making the necessary tests, no pump manufacturer will have this infor- 
mation unless it is required and paid for by the purchaser. The lack of 
this information means that in most cases when investigating the tran- 
sient hydraulic conditions in a pump discharge system, the pump char- 
acteristics in these regions must be estimated from the known curves 
of the very few pumps that have been run through the necessary tests. 

It would be of interest to have the author’s opinion on whether the error 
induced in making such assumptions could be of sufficient consequence 
to require the necessary tests be made. 

It is interesting to compare the transient hydraulic conditions in the 
Tracy pumping system as designed, with the system without any valves. 
Neglecting friction, which in this case will not induce appreciable error, 


and using the values of 2p, K; and 2h given by the author, we find 


that the pressure rise at the pump would have been I.40Hp or an in- 
crease of approximately 20 percent over the pressure rise obtained 
with the valve. This indicates how effective a valve of this type can be 
when properly operated. The pressure drop at the pump will be the 
same as with the valve, and will occur in 3.5 seconds after power fail- 
ure to the pump motors, and a complete reversal of flow in the pump 
discharge line will occur in about 10 seconds. This indicates that 
valves with slow closure time can be of no value in the control of tran- 
sient hydraulic conditions as the maximum pressure rise, the maximum 
drop and the reversal of flow in the pump discharge line occur in a very 
short time. Further comments by the author on his experience in the 
proper valve closure time would be of great interest to this writer. 

Quite frequently it is desirable to know the maximum reverse speed 
of a pump and to have a complete time history of the transient hydraulic 
conditions in a pump system after failure of power to the pump motors. 
If the author has such information in the form of curves similar to 
those shown in Figure 14, it would be a valuable contribution to this 
paper. 


| 

| 
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P. LINTON.!—The author has performed a most valuable service in 
reporting these full-scale tests on pressure surges after power failure, 
information which is unfortunately published only in very few cases. 
They show up clearly the value of a full analysis of transient conditions 
in leading to savings in surge control devices such as surge tanks and 
control valves. 

The time history of the pressure surge, Fig. 3, is most interesting 
and brings out clearly the results of the pressure surge diagrams, 
Figs. 11 and 12. Comparison between the measured and computed 
curve shows good agreement as far as the general shape is concerned; 
the lowest pressure reached is some 17 ft higher than the calculated 
one, while the subsequent pressure rise is 26 ft less than that computed. 
In both cases the measured pressures are more favourable than the 
calculated ones and the present tests therefore confirm that a careful 
analysis of the pressure surges leads to a somewhat conservative re- 
sult. This difference between calculated and measured results might 
be valuable as an additional margin of safety if it could be relied on in 
all cases, and it is important therefore to investigate some possible 
reasons for it. 

It can of course be taken that the instruments used were capable of 
following accurately the fast pressure oscillations which were expected; 
no details are given in the text and it would be interesting if the author 
could indicate briefly the type of instrument used and its probable ac- 
curacy under the actual conditions. 

The theoretical minimum pressure depends to some extent on the 
assumptions implicit in the surge diagram, Fig. 11. As noted in the 
text, the actual pipeline, which consists of a number of different pipes 
in series, has been replaced by the equivalent uniform pipeline deter- 
mined by equs. (13) and (14). This may have contributed to the low 
value of the minimum pressure since in a somplete diagram the line 
Ao A3.5 (Fig. 11) would be replaced by a series of lines staggered 


either to the right or to the left. It would be interesting to know 
whether the author has attempted to draw a complete pressure surge 
diagram to include all or at least most of the changes in the pipeline, 
and whether that has given a value closer to the measured one. 

Fig. 11 also indicates that the position of point Ags and with it the 


minimum pump pressure are determined largely by the shape of the 
pump head-discharge characteristic o@ in this region. A comparison 
with Fig. 5 shows that this represents a very considerable extrapola- 
tion beyond the design point, well past the end of the pump character- 
istic given in Fig. 4. It would be interesting to know the method by 
which the author has carried out this extrapolation, and whether the 
curves in Fig. 7a obtained by Prof. R. T. Knapp were used in this con- 
nection too. The writer would also like to know whether in the author’s 
opinion this part of the pump characteristic is likely to vary from one 
machine to another, particularly when the pump does not cavitate as in 
the present problem. 


1. A.M.I.E.E., A.M.LMech.E., Mem. A.S.M.E. British Hydromechanics 
Research Association, Harlow, England. 
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Additional and highly important results appear to have been obtained 
by the author in another test at Tracy pumping plant when only one pump 
of a pair discharging in a common pipe line was stopped. To the 
writer’s knowledge, no such experiments have been reported on an 
installation of this type, in which reverse flow through the pump occurs 
during slowing down before the butterfly valve has been closed. It 
would be valuable to have a comparison between calculated and meas- 
ured pressures for this case also, if this has been investigated by the 
author, together with some indication of the method by which the distri- 


bution of flow between the running pump and the one slowing down has 
been determined. 


| 
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DISCUSSION OF GRAPHICAL AND THEORETICAL 
ANALYSIS OF STEP-DRAWDOWN TEST OF ARTESIAN WELL 
PROCEEDINGS-SEPARATE NO. 362 


RAPHAEL G. KAZMANN,' A.M. ASCE.—This paper is a lucid and 
thorough treatment of well-hydraulics from an eminently practical 
standpoint. It will be of immediate value to all ground-water engineers 
who have occasion to design wells, prepare specifications for alternate 
designs of wells, or compare well-fields with other methods of water 
collection such as infiltration galleries. Furthermore it spells out in 
considerable detail our need for research to supply deficiencies in 
design-data. 

The total drawdown in a well is a major factor in determining the 
cost of bringing water to the surface of the ground. Well-life is a 
second factor of importance inasmuch as all calculations of annual 
cost are based upon this figure. The author has provided a basis of 
well-design which will tend to decrease drawdown and increase well- 
life and thus reduce the cost of water. 

Clear indication is given that wells should be designed as laminar- 
flow devices. Apparently, for any given well, as a certain critical flow 
rate is exceeded, the drawdown is affected disproportionately by factors 
which the designer can modify: well-diameter and screen area, for 
example. The significant flow rate appears to be Q,. 

The author has evaluated Q, for a screen 11-inches in diameter and 
30 feet long. This value, 0.75 cfs, is of great practical importance. It 
enables computation of the velocity of water through the screen, or 
through the formation immediately outside the screen. 

Although the total area of screen-opening is not given (and this 
would have been desirable), the porosity of material next to the 
developed well would probably range from 40 to 25 percent. Thus, the 
velocity of water just before it reached the screen surface (neglecting 
screen thickness) would range from about 0.022 ft/sec to 0.035 ft/sec, 
values significantly lower than figures generally cited in the literature. 

In this connection it might be noted that an increase in the area of 
screen opening beyond 40 percent of the total screen area is not likely 
to produce a desirable decrease in velocity, all other factors remaining 
the same. Ideally the percentage of total screen area in openings 
should be the same as, or slightly greater than, the porosity of mate- 
rial immediately surrounding the screen. A highly uniform material 
might justify a larger percentage of openings than a less uniform mate- 
rial. Percentage of opening would normally range from about 15 per- 
cent to about 35 percent of the total screen area. 


1. Consulting Engineer, Stuttgart, Arkansas. 
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It would be of value to the profession if the author could make avail- 
able a table listing the following information for the several tests 
mentioned in the paper, and especially those listed on page 362-6: Q., 


BD, C, P, F, &, Ty length of screen, and percent of total screen area in 


openings. Such a table would show “order of magnitude” figures for 
the upper limit of velocity of laminar flow towards wells. An additional 
column in the proposed table might be labelled “v _" the critical 


velocity below which laminar flow prevails from the aquifer to the in- 
ner surface of the well screen. 

The magnitude of v,, will probably prove of importance in designing 
for increase of well-life. Experience has shown that, in incrustating 
waters, the more intensively a well is pumped, the shorter its life. 
When a ground water contains calcium bicarbonate it tends to deposit 
lime on the screen and in the interstices of the material immediately 
surrounding the screen. It has been observed that, as this process 
continues, the effective area through which water can flow is reduced, 
the drawdown increases, the rate of incrustation is increased dispro- 
portionately, and the well must be treated to restore yield or must be 
abandoned. Present theory indicates that where the cone of depression 
is steepest, that is, where the reduction in the internal pressure of the 
water is most rapid, the deposition of calcium carbonate occurs most 
rapidly. Equation of the reaction is approximately 


Ca(HCO,) ----» CaCOg}+ CO,} + H,O 
3’2 3 2 2 


It follows, therefore, that elimination of turbulent flow in the vicinity 
of a well should eliminate the abrupt pressure drop associated with 
turbulent flow and thus retard this chemical reaction. Well-life would 
thus be increased and maintenance costs greatly lowered. 

Manifestly there are great opportunities for constructive work in the 
field of well-hydraulics under field conditions. Such studies should in- 
clude a series of tests in different places under differing hydrogeologic 
circumstances to yield conclusive observations on factors of immediate 
practical importance such as Q, and Vor In conjunction with such 


studies it would be desirable to carry on studies of water quality to 
determine, if possible, in what manner the rate of CaCO, deposition 
(and deposition of other minerals) is related to the steepness of cone of 
depression. Such studies would result in immediate and tangible bene- 
fits to well-users and aid in the conservation of power and critical 
materials. 

ABDEL-AZIZ I. KASHEF,' A.M. ASCE.—Mr. Rorabaugh raised in 
his paper a matter of major concern, namely, the indication “that the 
well radius maybe more important in well design than has been con- 
sidered in the past.” % 

Under the steady-state condition the discharge varies as log = and 


n 
this fact leads to the conclusion that the effect of well radius ¥. is not 


1. Lecturer, Ibrahim University, Cairo (Egypt), Formerly Research 
Fellow, Purdue University, U.S.A. 
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r, 
important because the term log = varies little as Th is changed. 
n 


Under the non-equilibrium condition, the well known formula arrived at 
by both Theis? and Jacob’ assumes that the pumped well has a theoreti- 
cally infinitesimal diameter or a hypothetical line sink. In order to 
apply this non-equilibrium formula to an actual problem, Jacob sug- 
gested the determination of the drawdown curve for a line sink placed 
at the center line of the actual well, then the intersection of the well 
surface with the drawdown curve determines the accurate drawdown at 
this surface. This method was based on the fact that the rate of with- 
drawl of water from storage within the well casing is negligible.* This 
means that the effect of well radius is considered negligible under the 
non-equilibrium flow towards wells. 

However, it had been shown by the writer® and others using numerical 
methods based on finite-difference equations that there is a certain ef- 
fect of the well radius that leads to a variation of the drawdown-—the 
discharge and other factors being constant--near the vicinity of the 
well casing compared to the Theis solution. 

An example had been worked out by the numerical method and the 
Theis equation, and both solutions conformed more or less precisely 
only beyond 50 feet from the well center. It was concluded that, “This 
result was expected since the Theis formula is derived on the assump- 
tion that the well is represented by a hypothetical line sink. This as- 
sumption results in infinite values at the center line of the well. The 
Theis curve will not intersect the well surface at a point that corres- 
ponds exactly to the actual drawdown. In order to do this the well cas- 
ing itself would have to be filled with the same kind of material as that 
through which the water percolates and not filled merely with water. 
The assumption of a line sink is not made in the numerical solution. 
The actual boundaries are considered.”* The data chosen for this 
solved example were the same as that applied in the Grand Island 
Tests®: discharge = 540 gallons per min., transmissibility = 8.866 sq. 
ft. per min., coefficient of storage = 0.217, radius of well = 1.0 foot, 
the initial piezometric head was assumed to be 100 ft. 


2. Theis, C.V., “The Relation Between the Lowering of the Piezometric 
Surface and the Rate and Duration of Discharge of A Well Using 
Ground-Water Storage,” Trans. Amer. Geophys. Union, Vol. 16, 
pp. 519 - 524, 1935. 

. Jacob, C.E., “On the Flow of Water in An Elastic Artesian Aquifer,” 
Trans. Amer. Geophys. Union, Vol. 21, pp. 574 - 586, 1940. 

. Jacob, C.E., “Flow of Ground Water,” Chapt. V in “Engineering 
Hydraulics,” edited by H. Rouse, John Wiley, 1950. 

. Abdel-Aziz I. Kashef, Y.S. Touloukian and R.E. Fadum, “Numerical 
Solutions of Steady-state and Transient Flow Problems—Artesian 
and Water-table Wells,” Engineering Experiment Station Bulletin 
No. 117, Purdue University, July 1952. 

. Wenzel, L.K., “Methods for Determining Permeability of Water- 
Bearing Materials,” U.S. Geol. Surv. W.S. Paper 887, 1942. 
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The maximum difference between both solutions—Theis and the 
numerical solution—was at the well surface. The drawdown in feet at 
this surface and after different time intervals are shown in the given 
table. It is apparent that the Theis solution gives greater drawdowns 
as compared to the numerical solution. The difference ranges from 
about 80% of the drawdowns calculated by the numerical method in the 
early stages of pumping to about 50% in the late stages. 

The writer’s opinion agrees to the conclusions arrived at by Mr. 
Rorabaugh concerning the importance of the radius on well design. 
Moreover, there is no doubt that under the non-equilibrium conditions 
the difference in drawdowns aforementioned is due to neglecting of the 
well radius in the derivation of the Theis formula. More investigations 
are in urgent need in this respect so as to clarify the matter and arrive 
at general conclusions regarding the effect of well radii. 


Drawdown (In Feet) at the Well Surface 


Time Elapsed 

Since Start of 

Pumping (in 

Hours) 

Drawdown 

Calculated by 6.742 7.640 7.903 | 8.090 
Theis Formula 


Drawdown 

Calculated by 

the Numerical 4.758 5.745 
Method 
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CORRECTIONS TO FREQUENCY OF EXCESSIVE 
RAINFALLS IN FLORIDA 
PROCEEDINGS-SEPARATE NO. 366 


Because additional data have been obtained for the rainfall at 
Jacksonville, Fla., Eq. 4e should be changed to 


231 


i 


This additional data will be included in the paper when it is printed in 
Transactions. 
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DISCUSSION OF A NEW CONCEPT OF FLOW 
IN ROUGH CONDUITS 
PROCEEDINGS-SEPARATE NO. 390 


V. L. STREETER,’ M. ASCE.—In attacking the difficult problem of 
classifying both artificial and natural roughness in terms of the resist- 
ance function and wall Reynolds number the author is to be strongly 
commended. The lack of adequate experimental results does not as yet 
permit complete evaluation of the new concepts. 

The author’s analysis of quasi-smooth flow makes the tacit assump- 
tion that the kinetic energy per unit time flowing past a section through 
a vortex is the energy consumption of the vortex. There is probably a 
very slow mixing of the flow of the main stream with the fluid contained 
within the vortices, as the rapid shedding of vorticity would result in 
much larger energy losses in the pipe flow. 

It is reasonable to assume that the rate of energy conversion to heat 
within a vortex in unit time depends upon the energy of the vortex. The 
expression for kinetic energy of a forced vortex ring is computed, then 
the vortex energy per unit weight of fluid flowing is developed. From 
this expression and the IIT experiments some conclusions as to energy 
dissipation within a ring vortex may be drawn. 

The kinetic energy contained in a ring vortex of circumference 7D 
and filament radius j/ 2 angular velocity w is 


16 


2 
1 1 P 
v2dm = fi 2mrdr = 12 pDCyv,,? 


0 


The number of vortex rings per unit length of pipe is 1/2j for square 
groove roughness with land thickness equal to groove width. The weight 
flow through the pipe per unit time is 


pevr D2/4. 


The kinetic energy in the vortex rings per unit weight of fluid flow- 
ing per unit length of pipe is 


cv 


As Vw/V should not vary substantially, the energy in the vortices per 
unit length per unit time is proportional to jV/(gD). 


1. Research Professor in Mechanics, Llinois Institute of Technology, 


Chicago. 
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Three sizes of square groove roughness were tested at IIT, with 
dimensions varying in the ratios 1:2:4. The test data, Fig. 7, indicates 
that the friction factor is practically the same for the three roughnesses, 
despite their large variation in relative roughness, and that they are 
constant for Reynolds numbers larger than 200,000. The energy in the 
vortices per unit weight of fluid flowing varies directly as the size j of 
the vortex and directly as the average flow velocity or wall velocity. 
Letting t be the fractional part of the energy converted into heat per 
unit time within the vortices, 


where k is a proportionality factor. Solving for f' the resistance coeffi- 
cient due to the vortices, 


jt 
Vv 

This shows that t must vary inversely as j and directly as V for constant 
f'. Hence the larger the vortex cross section (with correspondingly 
lower angular velocity) the less the rate of energy dissipation per unit 
of energy available in the vortex. Similarly the proportion of energy 
dissipated per unit of energy in the vortex increases directly as the 
wall velocity. 

WALTER RAND,’ A.M. ASCE.—The remarkable attempt, made in 
this paper, to analyze the flow over rough surfaces by taking into con- 
sideration the roughness characteristics h, A, and s demonstrates at 
the same time the complexity of the problem and the difficulty of de- 
riving equations that reflect all the aspects of the flow. The writer 
analyzed some of the equations related to the wake-interference flow in 
an attempt to find additional ways of interpretation and to reveal the 
significance of the numerical coefficients involved. 

The logarithmic velocity distribution, adopted by the author as the 
basis of analysis (Eqs. 3 and 4), can be expressed in the simple form 


a = 1 log, Zz 
v* k y' 


where y’ is the value of y at which the logarithmic curve intersects 
the y-axis. In the hydraulics handbooks y' is given by 


(26) 


(27) 
for the smooth pipe flow, and by 
(28) 


for the rough pipe flow, with m = 30 for Nikuradse’s sand roughness. 
Introduction of these values of y' into Eq. (1) will give the equations 


of von Karman for * and of (see Eqs. 4 and 5 in this paper). 


1. Asst. Prof of Fluid Mechanics, New York Univ., New York, N. Y. 
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Analogous equations, involving the longitudinal spacing A, can be de- 
rived for the rough flow by assuming that y' is proportional to A : 


A = ny' (29) 
An alternative form of Eq. 7 would then be 


1 1 
= s Be + (30) 


A 5 Be n (31) 


The author’s assumption that A is a universal constant would mean in 
this case that the ratio A/y' should be a constant for all types of rough- 
ness (provided k is constant). However, as the roughness height h 
alone seems not to be sufficient to characterize the roughness fully, it 
is also difficult to believe that A alone could be sufficient. 

In fact, the author gives the values of A ranging from 6.8 to 9.3. 
The corresponding values of n will then range from 15.2 to 41.6. Such 
a scatter does not well warrant the assumption of a constant A. It 
seems that A is rather a function of all the roughness characteristics 
(h, A, Ss, cross-sectional form of roughness elements) and could be 
constant for geometrically similar roughness. That means that the 
coefficient “1.75” in Eqs. 15 and 16 does not have a general validity 
and the value of resistance function for fully rough flow (high R,) 
should rather be written in the general form 


1 


. 1 1 To 1 1 
f d 


for two-dimensional flow. 

The author assumes that k is a constant (k = 0.4) for the total cross- Fs 
section of the flow in the case of fully rough flow as c reaches zero at 
high values of the wall Reynold’s number. In hydraulics the coefficient 
k is also accepted to be 0.4 in the case of smooth flow with the laminar 
sublayer covering the roughness elements. If k would be the same 
also for the transition zone between the smooth and fully rough flow, 
the conclusion would be that a variable n in Eqs. (32) and (33) could 
give the law for the resistance function in this zone. However, relying 
upon the results of velocity measurements (see Fig. 1) the author adopts 
a variable k (by 1/k = W) for the wall zone (Eq. 8), to derive the 
Eqs. 11 and 12. These equations reflect the correct trend of the 
resistance function for the transition zone in the case of uniform 
roughness. Nevertheless, even as the logarithmic velocity distribution 
curves (Fig. 1) show a tendency for a more pronounced slope near to 
the boundary, this tendency is artificially stressed by measuring y 
arbitrarily from a datum through the crests of the roughness elements. 
The choice of this datum would be inconsistent with the assumption of 
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the author that at high values of Ry, the coefficient k is valid for the 


total velocity profile (c = 0). This assumption requires that the veloc- 
ity should be zero at a finite distance y' over the crest of roughness 
that does not correspond to reality. In fact, various investigators 
(Nikuradse, Schlichting, Einstein, and others) have found that for fully 
rough flow a practically straight line velocity distribution can be 
plotted on logarithmic paper choosing the datum considerably below 
the crest of the roughness elements, rather somewhere near the 
centerline of these. Of course, the values of c and W will depend on 
the choice of the datum from where y has been measured. Considering 
also the fact that the actual k can be different from 0.4 (in Fig. 1 the 
values of k can be found to be 0.33; 0.42; 0.35), the use of the Eqs. 11, 
12, 15 and 16 seems rather difficult. 

The author uses the Eqs. 15 and 16 to derive the limiting values of 
resistance function (3.52 and 2.64 on page 390-11). However, for the 
large value of c in cA = ry the Eqs. 11 and 12 should be used. The 
—_ of resistance function for Ww = 0 and a. = 1 would then be 
ve for both types of flow (pipe and a. BEN For a value 


A = 8.5 the value of the resistance function would be 3. However, the 


limit we , with a value of A determined for fully rough flow can 


have no significance because a velocity distribution with W = 0 over 
the total cross-section is not physically possible. The flow with the 
maximum value of the resistance function is very near to the smooth 
pipe flow, that will define the limit of the resistance function, and 

W = 1/k and c = 0 would be a more realistic assumption, even if it is 
not consistent with Eqs. 11 and 12. Moreover, the transition from 
smooth flow to fully rough flow depends possibly as much on the ratio 


a (where 6' is the thickness of the laminar sublayer) as on the ratio 


A , 
5! and a wall Reynold’s number r/h could be more relevant than 


h 
RVE 


Toh to determine the variation of the resistance function: in the 


transition zone. 

The writer hopes that these considerations are a farther illustration 
of the difficulties encountered in the study of flow over rough surfaces. 

HARRY H. AMBROSE,! A.M. ASCE.—A review of data, published in 
the author’s references for small commercial and artificially- 
roughened pipes and elsewhere (1) for large conduits, supplies convinc- 
ing evidence to the fact that the effect of surface roughness character- 
istics upon the resistance to turbulent flow in conduits is an exceedingly 
complex phenomenon. The author is to be commended for his presenta- 
tion of a concept of this phenomenon that helps to clarify many of its 
peculiarities—particularly that of the rising resistance coefficient 


1. Associate Prof., Civ. Eng. Dept., Univ. of Tennessee, Knoxville, 
Tenn. 
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curves associated with “wake-interference” flow. The characteristics 
of the three basic types of flow are well defined and appear to satisfy 
the trends of the available data. 

The author has concluded that wake-interference flow, with its ten- 
dency toward constant resistance coefficients for moderate to high 
Reynolds numbers, is experienced infrequently in commercial pipes 
and that isolated-roughness or skimming flow is more likely to occur 
with attendant dependence of the resistance coefficient upon the viscous 
characteristics of flow for moderate, or even high Reynolds numbers. 
The writer anticipates a vigorous reaction to these conclusions but 
there is to be found much substantiating evidence. The appearance in 
the literature of data showing resistance curves closely paralleling the 
smooth-pipe curve has, in the past, been rather summarily dismissed 
and the surfaces yielding such resistance curves have been character- 
ized as “wavy” This description appears hardly adequate in view of 
the similarity of results obtained from such widely varying surfaces 
as towing planes coated with ship-bottom paints (2), a large steel pen- 
stocks (3) and nominally smooth pipes with artificial depression-type 
roughnesses. It is apparent that more serious consideration must be 
given to isolated-roughness and skimming flows. 

Although the writer finds phenomenological satisfaction in the three 
basic types of flow he disagrees with certain details of the relationships 
developed by the author. Every investigator of flow in rough pipes has 
been faced with the problem of deciding upon the correct diameter to 
be used in calculating the resistance coefficient and Reynolds number, 
and as a reference for velocity distributions. The choice is very im- 
portant both in relating velocity distributions at wall distances within 
the magnitude of the roughness dimensions and in determining the 
resistance coefficient. In the latter usage, the value of the resistance 
coefficient varies directly as the fifth power of the diameter. The 
author has chosen to measure the diameter between the crests of the 
roughness elements at the opposite walls. This choice appears logical 
for closely-spaced elements. In the opinion of the writer, however, the 
correct diameter depends upon the type of roughness—particularly, 
upon whether flow in the space between the roughness elements can be 
considered as part of the mean flow. As an extreme, the author’s rule 
would most certainly be incorrectly applied in large-scale isolated 
roughness consisting of protruding gaskets at flanged joints. The 
choice of diameter is also of particular significance in testing whether 
resistance effects of non-uniform elements are additive. 

To be truly general, the resistance equations must allow for transi- 
tion from smooth-pipe flow to each of the three basic types of rough- 
pipe flow and also for transition, whether abrupt or gradual, between 
the three basic rough-pipe flows. The author has not considered transi- 
tion from smooth-pipe to isolated- roughness or skimming flow although 
this is at least as important as the wake-interference transition. Tests 
by the writer upon roughness elements consisting of cylindrical depres- 
sions in a surface that is otherwise smooth have yielded well-defined 
transitions from smooth-pipe flow to skimming flow. Three such 
transitions are indicated in Fig. 1. Noteworthy is the appearance of 
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the transition of curve B which gives every indication of the approach 
to constant resistance coefficient to often assumed for rough pipes. 
The marked downward trend of f for higher values of R, however, is 
very closely parallel to the smooth-pipe curve and thus reveals the 
transitional nature of the preceding portion of the curve. (The symbols 
used herein are identicai with those defined by the author unless other- 
wise noted). 

The writer is at a loss to understand the author’s statement that for 
isolated-roughness and skimming flows the value of f is dependent upon 
R but not upon rp. Perhaps independence from the parameter To/ h is 


meant; however, the heretofore unpublished data of Fig. 1 show that 
To/ h is a controlling factor in the transition range. Values for the 


roughness-element geometry for the three curves of Fig. 1 are given 
in Table I. 


Table I 
Geometry of Roughness Elements 
Curve T/A A/h d/h ay/A, 
A 7.75 31 4 1 0.047 
B 7.75 31 4 2 0.189 
Cc 3.75 15 4 2 0.181 


(peripheral and longitudinal spacings are equal) 


(The symbol d represents the diameter of the depression, a,/A, is the 


net ratio of depressed area to total surface area, and the other symbols 
are as defined by the author.) 
The value of f cannot be independent of Tg inasmuch as comparison 


of curves B and C shows the difference in the transition portions to be 
due to To/ h (or Tp/¥)- Comparison of curves A and B indicates dis- 


agreement with the author’s Eq. (23) since the additive f-value appears 
to be closely proportional to ag/A, and independent of A/h. The claim 


may be made that the discrepancy is due to a variation in the product 
CwY yy» but this would certainly appear fortuitous in view of the agree- 
ment of the above-mentioned trends with the results of tests on fourteen 
other surfaces (which, for lack of time and space, are not included). 

For wake-interference flow, two distinct regions of velocity distribu- 
tion are defined—the central or core region and the wall region. The 
author has assumed velocity distributions for each region in the develop- 
ment of the resistance equations. Certain discrepancies are to be noted. 
Inasmuch as Eq. (14) is independent of the Reynolds number of the 
mean flow it can be valid only for the range proposed by Nikuradse— 
that of “completely-turbulent”flow. Likewise, it is not obvious that A 
is a constant and independent of both R and Ry unless the resistance 


equation follows the quadratic law. In spite of these facts, the distribu- 
tion of Eq. (14) is assumed to apply to the core region over the entire 
transition from smooth-pipe to rough-pipe flow. 
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If, as asserted by the author, the value of A is constant the core 
distribution of velocity will be governed only by the relative roughness 
spacing, whereas the wall distribution will be a function of both the 
type of roughness and the wall Reynolds number. For the assumed 
velocity distributions to apply over the entire transition range, there- 
fore, the wall region must contain all of the viscous effects. Since for 
a smooth pipe the entire velocity distribution is dominated by the 
Reynolds number, the following confusing picture presents itself. 

At the smooth-pipe end of the transition range the wall region en- 
compasses the entire turbulent flow, the velocity distribution being 
governed by the wall Reynolds number and independent of the rough- 
ness characteristics. As the Reynolds number is increased there ap- 
pears a core distribution at the center of the pipe governed solely by 
the relative roughness spacing although the wall velocity distribution 
enclosing the roughness elements is still primarily a function of the 
wall Reynolds number. Further increase of the Reynolds number 
causes the wall region to shrink and at the same time to be governed 
increasingly by the roughness characteristics rather than by the viscous 
effects. Finally, the wall region disappears and completely-turbulent 
flow ensues. In brief, if A is constant the effect of the roughness ele- 
ments must begin at the pipe center and radiate toward the wall as the 
Reynolds number increases. 

It does appear reasonable, however, that the velocity distribution 
should be divided into distinct zones which indicate the relative extent 
of influence upon the distribution of each governing quantity. As a mat- 
ter of fact, Ross (4) has shown this to be necessary for the more simple 
case of turbulent flow over smooth surfaces. Ross divides the entire 
turbulent zone into an “inner turbulent” region which is said to be 
governed by a universal “law of the wall,” an “outer turbulent” region 
governed by the spatial history of the turbulence, and a “blending” 
region which separates the inner and outer regions. Analogically, to 
the writer a more satisfying concept of the velocity distribution for 
transitional flow would be: (1) smooth-pipe velocity distribution as 
long as a laminar sub-layer persists; (2) with sufficient increase in 
the Reynolds number, generation by the roughness elements of eddies 
which initially would be weak and whose influence would be confined to 
the wall or inner turbulent region; and (3) with further increase in the 
Reynolds number, an increase in strength and extent of the eddies until, 
finally, the whole turbulent velocity distribution is dominated by the 
roughness elements. It is asserted that the influence of the surface 
roughness must extend from the boundary toward the center rather 
than the converse. 

The foregoing criticisms do not lessen the writer’s appreciation of 
the general aspects of the author’s concept of turbulent flow in rough 
pipes—namely, that three basic types of flow are possible and that the 
spacing of the roughness elements is the controlling factor as to which 
type of flow will attain in a given case. It is believed that the profes- 
sion would do well to reanalyze existing data in the light of this concept 
and to use it as a guide in future research on the subject. 
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HENRY M. MORRIS,' A.M. ASCE.—It has been somewhat surprising 
to the writer that, in view of the rather sharp departure of the new con- 
cepts from previous methods in a field of wide interest in hydraulics, 
there has been such relatively little criticism of the paper. It is hoped 
that this fact indicates the new concept will prove generally acceptable 
and useful to the profession. The comments and suggestions of the 
discussions by Messrs. Streeter, Rand, and Ambrose, together with the 
generally favorable attitude taken toward the paper, are greatly ap- 
preciated by the writer. 

Dr. Streeter suggests an alternate analysis of quasi-smooth flow 
friction factors, which indicates that the proportion of groove vortex 
energy actually being converted into heat varies inversely as the groove 


width and directly as the flow velocity, assuming constant friction factor. 


However, it appears to the writer that there is a dimensional incon- 
sistency in this analysis. Dr. Streeter derives an equation for the total 
kinetic energy, in ft.-lbs., contained in the groove vortex, and then 
divides the result by the weight-rate of flow, in lbs. per sec. The 
resultant expression, which Dr. Streeter calls the kinetic energy in the 
vortex rings per unit weight of fluid flowing per unit length of pipe, is 
thus actually the vortex energy. for each lb. per second flowing, for a 
unit length of pipe. The dimensions for his expression are therefore 
(ft.-secs.) per ft., and the quantity seems to have no pertinent physical 
significance as it stands. Also, it seems questionable to attribute the 
entire flow friction factor, as measured in the LI.T. tests, to the energy 
of the vortex. A substantial part of the energy consumption of the flow 
is due to the normal “smooth-pipe” action, whereby vorticity is gener- 
ated and shed into the main flow at the pseudo-wall consisting of the 
land areas of the threads and the upper limbs of the groove vortices. 

It is only the additional energy loss, above the normal smooth-pipe 
loss, that is attributable to the groove vortices. 

It seems reasonable to assume, as the writer did in the derivation 
of Eq. (22), that the entire vortex energy per unit time is being con- 
tinuously converted into heat, under equilibrium conditions. However, 
the vortex itself is not dissipated, as it continually receives an equiva- 
lent supply of power through attrition of the bulk flow with its upper 
limb. 


Lafayette, La. 
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In reference to Eq. (22), it should be noted that there are two typo- 
graphical errors. Thus, the term R should be Ry the hydraulic 


radius, and the term V should be y, the specific weight. 

Professor Rand criticized the theory of wake-interference flow but 
the writer believes that the theory in its essentials is adequate to 
describe this type of flow, especially since it appears to be substanti- 
ated in general by the available experimental data. The phenomena 
are very complex and the theory and equations as derived involve a 
number of simplifying assumptions, but they nevertheless appear more 
adequate to explain the phenomena involved than do previous methods. 

The roughness element spacing seems more reasonable as a physi- 
cal measure of the roughness than does the roughness height, especially 
in wake-interference flow, in which the turbulence phenomena originate 
at or near the crests of the elements and thus can be affected very 
little by the height of the elements. This also is sufficient argument 
for taking the datum from which radial distances are measured at the 
crests of the elements, rather than at some altogether arbitrary datum 
selected purely to make the same logarithmic velocity distribution 
slope apply throughout the pipe. Both physical considerations and ex- 
perimental measurements have shown that these slopes are not the 
same in the wall and central zones, but that there must be a greater 
intensity of turbulent mixing than is reflected by the von Karman k in 
the wall zone than in the central zone. Even the fact that this procedure 
must ignore a small transition zone near the wall does not make it as 
physically unrealistic as would the arbitrary extension of the central 
zone velocity distribution equation not only through the wall zone of 
abnormal turbulence but even into the dead water regions between the 
elements, as is true of the methods defended by Professor Rand. Added 
to these considerations is the fact that the new concept and procedure 
provides at least a partial explanation for the peculiar characteristics 
of the friction factor and resistance function curves for this type of 
flow, which the older methods did not, and provides a very effective 
basis for the correlation of data from many types of roughnesses, as 
shown in the paper. 

Professor Rand also criticizes the assumptions made by the writer 
in determining the maximum possible value of the resistance function. 
It is recognized, of course, that these assumptions are extreme and 
actually represent physical impossibilities. However, they are quite 
reasonable as limiting conditions, and actually come close to being 
realized for low wall Reynolds Numbers. Both measurement and 
physical considerations indicate that YW approaches 1/k and c approaches 
zero at high values of Ry» rather than at the point of initiation of wake- 
interference flow. 

The question of the location of the datum for calculating radial dimen- 
sions is also discussed by Professor Ambrose. He rightly emphasizes 
the importance of this question, as it has a very great effect on the fric- 
tion factor and velocity distribution. It seems evident that the rough- 
ness crests provide the best datum in the cases of skimming flow and 
wake interference flow. The remarkable correlations of experimental 
data from many different sources and on many types of roughnesses, as 
discussed in the paper, are possible only when the datum is chosen thus. 
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However, Professor Ambrose is probably correct in pointing out that 
isolated roughness flow, particularly when the roughness elements are 
widely separated, and when they have large radial heights, should be 
described in terms of the wall-to-wall diameter. The case of the 
transition region between isolated roughness and wake-interference 
flow, and the associated problem of the correct datum, is one which 
needs further investigation. 

The new data on skimming flow obtained by Professor Ambrose are 
very interesting and provide a welcome extension to the I.I.T. and 
other data on this type of flow. An item of interest is that the f-R 
curves for the LLT. roughnesses, which were strip-type roughness 
elements, separated from the smooth-pipe curve at much lower values 
of Reynolds Number than did the curves for the Ambrose spot-type . 
roughness elements. The inference is that, in skimming flow, the 
larger the ratio of wall area to depression area, the larger will be the 
value of R at which the depression vortex begins to contribute substan- 
tially to the friction factor. 

The writer’s statement that, in isolated-roughness and skimming 
flow, the friction factor was independent of the radius was not, of 
course, meant to be taken as excluding the effect through the Reynolds 
Number, which was explicitly stated to be functionally related to f. 

It was meant to say that the friction factor was independent of the rela- 
tive roughness and the relative roughness spacing for these types of 
flow. Apparent variations of friction factor with relative roughness 
as obtained on commercial pipe tests can probably be attributed to in- 
correct choice of radius for computation of friction factor. This had, 
in fact, previously been suggested by Piggott,’ who, in his extensive 
studies of commercial pipe friction factor data, noted that the f-R 
curves for a given roughness type could usually be made to coincide, 
even for pipes of different radii, provided the proper datum was 
chosen. Except for the unusual conditions mentioned previously, this 
would most likely be the datum through the crests of the roughness 
elements. 

The cylindrical depression spot-type roughness elements tested by 
Professor Ambrose are different from any coming previously to the 
writer’s attention, and it is difficult to classify them. His curve A, 
for which the ratio d/h (corresponding to the writer’s j/h) was 1,seems 
certainly to follow the trends indicated in Eq. (23) for skimming flow. 
However, this equation was derived primarily for a groove-type depres- 
sion extending around the pipe periphery, and reflects the energy expen- 
diture in maintaining a stable ring vortex. This is a somewhat different 
situation from the vortices that would be formed in cylindrical spot 
depressions. Furthermore, the roughnesses represented in his curves 
B and C had d/h values of 2 and, as indicated by the writer, it is 
doubtful whether the depression vortices are maintained as stable 
entities when the groove width becomes materially greater than its 
depth. 


1. R. J. S. Piggott, “The Flow of Fluids in Closed Conduits,” 
Mechanical Engineering, Vol. 55, 1933. 
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Consequently it is doubtful that skimming flow was produced at all 
on roughnesses B and C. The typical form drag phenomenon of a 
continued generation of vortices in the wake would certainly extract 
more flow energy than the stable vortices of skimming flow, and this 
would most probably account for.the increase in friction factor as 
between roughnesses A and B, rather than the greater proportion of 
depression area to wall area, as suggested by Professor Ambrose. 

It is therefore likely that, whereas roughness A produced skimming 
flow, roughnesses B and C produced either wake-interference or 
isolated roughness flow, with wakes being produced at both extremities 
of each depression. In view of the peculiar character exhibited by 
Curve B, there is the possibility that roughness B_ tended to produce 
wake-interference flow at low Reynolds numbers, changing to isolated 
roughness flow at higher Reynolds numbers, with the result that 
curves B and C essentially coincide thereafter as would be indicated by 
Eq. (19) for isolated roughness flow. 

Professor Ambrose mentioned other experimental data which he had 
obtained. It is hoped that, when these are published, they will be able 
to clear up these problems. 

The writer believes Professor Ambrose, in his remarks about the 
velocity distribution equations in wake-interference flow, may not have 
grasped the full significance of the wall region of abnormal turbulence. 
It is not strictly correct to speak of this sort of flow as transitional 
between smooth-pipe and rough-pipe flow, or of the latter as being 
characterized by “complete turbulence,” because such terminology 
tends to imply that this transitional flow is not itself fully turbulent, 
but partly controlled by boundary layer effects remnant from smooth- 
pipe action. Actually, the energy losses in this type of flow are propor- 
tional to a power of the velocity greater than 2, as shown by the rising 
f-R curve, and thus the turbulence in the wall region should be described 
as abnormal turbulence, as distinct from the normal turbulence charac- 
terized by the quadratic resistance law, and also by the von Karman 
constant k. There is, of course, a short transition from smooth-pipe 
phenomena to these phenomena of abnormal turbulence; this transition 
is represented by the portion of the “resistance function—wall Reynolds 
number” curve which gradually diverges from the smooth-pipe curve, 
finally becoming horizontal. The long rise in the curve which then 
follows, finally becoming horizontal again, is the zone of abnormal 
wall turbulence. 

The wake-interference flow theory and equations, as developed in 
the paper, deal of course only with the part of these phenomena repre- 
sented by the rising and ultimately horizontal f-R curve. Isolated 
roughness flow would occur at lower Reynolds Numbers, represented 
by the portion of the transition curve from the point of departure from 
the smooth-pipe curve to the low point in the dip; when wake-interference 
flow is initiated, the f-R curve begins to rise. This reflects the fact that 
the wall zone has been transformed from one controlled primarily by 
phenomena in the laminar boundary layer, which by then has vanished, 
to one of abnormal turbulence caused by the “wake-interference” 
phenomenon. The velocity distribution near the wall in smooth-pipe 
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flow is characterized by a slope corresponding to a number greater 
than 1/k; the wall velocity distribution in wake-interference flow has 
a slope W which is less than 1/k, caused by the abnormal intensity of 
turbulence in that region. 

It has been found and generally accepted that the dimensionless 
velocity distribution equations in the central region of the pipe have the 
same slope (corresponding to the von Karman k) regardless of Reynolds 
Number and for both smooth and rough walls. This is the zone denoted 
by Ross (as cited by Professor Ambrose) to be the outer turbulent zone, 
governed by the spatial history of the turbulence. In this zone, the 
vorticity generated at the wall or elsewhere, by whatever means, 
whether at a laminar boundary film or by wake vorticity, has been 
broken down by mutual attrition and momentum transfer to “normal” 
turbulence. It is surely reasonable that the central zone turbulence 
would also be “normal” during the transition throughout the full range 
of wake-interference flow. 

As soon as the transition from smooth flow to the quadratic law 
flow is complete (indicated by the low point of the dip on the f-R curves), 
then true wake-interference flow begins. From this point on, the wall 
zone turbulence is governed by the roughness elements entirely (more 
precisely by the element spacing and the form of the wakes behind 
them, which latter is a function of the geometric form of the elements 
throughout the transition). It seems reasonable that, for this regime, 
the term A should be a constant as well as k. Reasons for this were 
discussed in the paper, as well as the experimental evidence that it 


was so. 
Thus, it appears to the writer that most of Professor Ambrose’s 


discussion of the phenomena of transition from smooth to rough pipe 
flow are quite correctly applicable to the first stage in the transition, 
but not to the second stage. At the initiation of this stage, the influence 
of the wall roughness has already extended completely to the center of 
the pipe, and the boundary constant A has become established as a 
characteristic of the core velocity distribution, although only in the 
central pipe zone has the wall turbulence been reduced to the normal 
turbulence characterized by k. Then, as the Reynolds Number increases 
further, the central zone extends toward the wall, both by shifting of the 
boundary of the two zones toward the wall and also by more rapid reduc- 
tion of the abnormal wall turbulence to normal turbulence even in the 
wall zone, as evidenced by the approach of the velocity distribution 
slope in the wall zone to 1/k, the slope of the core velocity distribution. 
There is no question but that the phenomena of flow in rough pipes 
are many and complex. There is much room for additional investigation, 
especially in the various boundary zones. The present writer hopes that 
the new concepts presented in this paper will help to clarify at least 
some of these phenomena and point the way to further research, and at 
the same time will prove useful in practical engineering design proce- 
dures. Again he wishes to thank the writers of the discussions for 
their interest and very material contributions to that end. 
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DISCUSSION OF BACKWATER EFFECTS OF OPEN 
CHANNEL CONSTRICTIONS 
PROCEEDINGS-SEPARATE NO. 413 


PAUL V. HODGES,' M. ASCE.—A very useful and simple method for 
the determination of backwater effects caused by different types of 
channel contractions is presented in this paper. The base curves and 
coefficients are defined in the authors’ reference (2). When the drop 
in water-surface, Ah, for a particular contraction has been determined 
by field observation, the actual backwater, h¥, caused by the contrac- 
tion, can then be determined from base curves similar to the authors’ 
Figure 9, which would give the backwater ratio h)/ Ah for different 
degrees of contraction. 

According to the tests made in the investigation, it is shown that the 
variation of L/b has no effect on the backwater ratio. In authors’ refer- 
ence (2) the coefficient C is shown to vary with L/b, and in the authors’ 
Figure 12 the backwater ratio is shown to vary with coefficient C. 
Consequently it appears that the backwater ratio should also vary with 
L/b. 

The base condition is defined as a vertical-faced constriction with 
square-edged abutments, and the curves in authors’ Figure 9 evidently 
refer to this condition only. It is inferred that the relation between 
h}/Ah and m would be different for each type of constriction. 

By assuming different values of Ah and y, for a particular type and 
size of bridge opening, the discharge Q can be computed by the equation 
in authors’ reference (2), and thus the relation between backwater h* 
and the discharge can be determined to fit the assumed conditions. 
Under actual conditions of flow, however, the relative values of Zh and 
and y, would be governed by channel characteristics, both upstream and 
downstream from the contraction, and consequently the discharge could 
be computed only after the actual determination of Ah and y,. 

This method, as presented by the authors, could often be used to 
compute the backwater from channel obstructions caused by bridge 
piers. Backwater formulas, based on many tests, have been derived 
for the computation of backwater effects of bridge piers, of which the 
Nagler formula appears to be the most satisfactory. See writer’s 
references (5) and (6). Using the authors’ notations, Nagler’s formula 
is 


Q = 8.02 Kybly, - 0V2/2g] [h* +B V3/2¢f 


The values of coefficients 6 and 8 are given in references (5) and (6). 
The relations between drop in water surface Ah, backwater h*, 
degree of contraction, m, and discharge as computed for a certain type 


1. Formerly with Corps of Engineers, Denver, Colo. Retired. 
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and size of bridge opening, are given in the following table. This table 
shows the computed value of K)> which when used in Nagler’s formula 


will give the backwater hj. 


Writer’s Table 1.—Relation Between Contraction, 
Backwater and Discharge 


Ah ht/Ah hi 


b = 90 feet n = .025 
L = 30 feet 

y, = 15 feet 

A, = 1350 sq. feet 


Q = Discharge as computed by the equation and coefficients 
given in authors’ reference (2). 


The authors indicate that the backwater ratio is not influenced by a 
variation in the paramters y,/} and L/b, but only by the degree ‘of 
contraction m. Under these conditions the method presented in this 
paper would have wide and easy application compared tu earlier 


methods which are much more complicated and involved with discharge 
formulas. 


References 


5) Woodward, S.M., and Posey, C.J., “Hydraulics of Steady Flow in 
Open Channels,” Chapter X, p. 127. John Wiley and Sons, Inc., 1941. 


6) Yarnell, D.J., Bridge Piers as Channel Obstructions, U.S. Department 
of Agriculture Technical Bulletin No. 422, 1934, 


CARL F. IZZARD,’ A.M. ASCE.—This paper,together with the com- 
panion paper by Kindsvater and Carter (1), reports the first comprehen- 
sive effort to obtain model experimental data from which to develop a 


1. Chf., Hydr. Research Branch, Bureau of Public Roads, U.S. Dept. of 
the Interior, Washington, D. C. 
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20 .34 .34 .886 12,900 
20 34 .68 .916 19,100 
30 45 .861 10,900 
30, 45 .90 .884 15,900 
40 .54 .840 9,470 
40 .54 1.08 .860 13,800 
50 .62 .62 8,510 
50 .62 1.24 .822 12,400 : 
Type.—Vertical face constriction with square-edged abutments. 
y, = y;+4h 
L/b = .333 
y,/b = .167 
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method for computing backwater caused by a channel constriction such 
as a bridge. The contribution reflects the judgment of the U. S. Geolo- 
gical Survey that small-scale hydraulic model tests offer a sound means 
of exploring the significance of the many variables involved in the prob- 
lem of indirect measurement of discharge by the contracted opening 
method. The research has been intelligently organized and aggressively 
conducted. The reports published to date (including this paper) have 
been confined to conclusions drawn from the experimental data. While 
not stated in the paper, the Survey is making field measurements of 
backwater and discharge as rapidly as possible in order to verify the 
results of the model tests. 

One of the outstandihg values of the work accomplished is the knowl- 
edge of where to take significant measurements of contracted flow. The 
verification of the laboratory results by prototype measurements will 
be followed with great interest not only by engineers who want to com- 
pute discharge for a known backwater but also by engineers who must 
estimate backwater for a proposed bridge or other constriction of an 
open channel. 

This discussion will be limited primarily to the problem from the 
viewpoint of the highway engineer since this is the field in which the 
writer is engaged. 


Viewpoint of the Highway Engineer 


The primary objective of the highway engineer is to build a structure 
which will safely handle the vehicular traffic at a minimum annual cost. 
That cost should include not only capital and maintenance cost but also 
costs for less tangible items such as flood risk (defined as flood damage 
cost multiplied by probability of occurrence). The flood damage cost is 
directly related to the amount of backwater both with respect to damage 
to abutting property, which sometimes entails litigation, and to destruc- 
tion of parts of the highway itself. The frequency of a flood which over- 
taxes the highway structure determines the period of time over which 
the probable damage must be charged off as flood risk cost. Depending 
on the nature of traffic and average daily count the flood risk cost may 
also include substantial charges for interruption to, or interference 
with, movement of traffic. Where approach roadways are designed for 
overflow, backwater and duration of flcod as well as frequency are 
important. 

One important distinction must be noted between the objectives of 
the hydrologic engineer and the highway designer. The former is ex- 
pected to achieve a fairly high standard of accuracy in his estimate of 
the flood discharge as computed from backwater and that answer is his 
end result. The highway designer, on the other hand, puts the computa- 
tion in reverse and wants to know how much backwater can be expected 
for floods of various frequencies whose magnitudes can probably be 
estimated not closer than twenty percent plus or minus (unless a gaging 
station with long record happens to exist somewhere nearby on the 
stream). Obviously then the highway engineer does not need to work to 
the close tolerances expected of the engineer who is gaging streams. 
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This distinction should be borne in mind to avoid creating the impres- 


sion that the writer is unduly simplifying a complex problem, or casting 


reflections on the U.S. Geological Survey for being too meticulous. At 
this stage of development the research must be as thorough as possible. 
The following analysis was possible only because of the cooperation 
of the authors who made available all the basic data recorded from the 
tests conducted at the A. N. Talbot Laboratory of the University of 
Illinois. 


Analysis of Data for Base Condition 


A constriction placed in an open channel requires that the velocity in 
the constricted section be greater than in the normal section. The 
forces producing this change in momentum neglecting forces to over- 
come normal resistance are the unbalanced forces represented by the 
total pressure at some upstream section less the total pressure against 
the embankment constricting the flow. Since the water elevations 
against the embankment are not known with any precision and are de- 
creasing where flow is accelerating, exact mathematical analysis in this 
reach is out of the question, although it does work fairly well in the de- 
celerating reach below the constriction as demonstrated in the paper. 

This concept does lead to an hypothesis that the increased depth in 
the approach section relative to the normal depth of flow should cor- 
relate with the change in momentum. The change in momentum ob- 
viously is dependent on the conveyance of the constricted section rela- 
tive to the conveyance of the approach section, which is essentially the 
same as (1 - m). A convenient form of plotting the experimental data 
as a means of testing this hypothesis is shown in Fig. 13 in which the 
ordinate is the depth ratio first mentioned. 

The abscissa in Fig. 13 is the velocity head in the constricted sec- 
tion relative to the “normal”depth. In this case, since the laboratory 
flume was level, “normal” depth means the depth Yp at the location of 


the constriction with the constrictions removed. The velocity, of 
course, is the discharge divided by the area at this “normal” depth. 
The difference in Y, at Sections 1 and 3 can be considered negligible 
in these tests. 

The experimental data for one size of vertical-face constriction with 
square-edged abutment (base condition) define three curves for three 
values of roughness. The maximum difference between extreme values 
of n is on the order of 2.5 percent of the normal depth, and is therefore 
not important from the highway viewpoint. Also the length of constric- 
tion L measured along the flow has little effect so long as L is less 
than one-third of the width of opening b. The curves also apply 
reasonably well to the non-rectangular section (Fig. 5) for a comparable 
value of m as shown by the two points so designated in Fig. 13. 

For other values of b the change in momentum for a given “normal” 
flow would be different. Figure 14 shows curves for values of m of 0.21, 
0.41, 0.61 and 0.81; no ordinate scale is indicated since the object is to 
show the relative position of the curves and not to provide a working 
chart. The curve for m = 0.41 corresponds to an intermediate rough- 
ness value in Fig. 13. 
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Application of Analysis 


The apparent correlations obtained from analysis of the data from 
the scale model are encouraging but not conclusive. The indications 
are that a graph such as Fig. 14 might be used as a simple means of 
estimating relative backwater caused by bridges of different length. 

A separate set of curves might be developed for each of the commonly- 
used bridge types. Adjustment factors for skew, eccentricity, and piers 
might be computed. 

Comparisons of backwater computed by this short-cut method and by 
the procedure recommended by the authors for certain hypothetical 
cases have given reasonable agreement. The agreement is better if the 
depth ratio at the upstream point is considered to apply to the mean 
depth at that section rather than the depth above the stream bed. 

The validity of any computational method based on the model tests 
cannot be firmly established until prototype measurements of backwater 
and discharge are available for verifications. It may be noted that 
actual stream crossings will rarely involve depths of flow as deep rela- 
tive to the flow width as used in the model test. 

While the use of a level flume simplified testing at both the labora- 
tories, the effect of stream slope on backwater is obscured. It seems 
reasonable to assume that the increment in head loss over the normal 
resistance loss in the approach reach with no constriction is included 
in the empirical evaluation of the ratio y;/ Y,: This is a weak point in 
the writer’s analysis since the model is obviously not geometrically 
similar to any usual prototype. The separation of losses as provided 
by the authors may be more nearly correct for some situations. 


The Backwater Ratio 


The backwater ratio as defined by the authors is apparently indepen- 
dent of the Froude number. The ratio y,/ Yin in Fig. 13 is a function of 


F for the constricted section at “normal”depth. This follows because 
the abscissa is 


/ = F2/2. 


Since y,+h* = (y,/y.) y.. by definition, it also is a function of F. 


Further if the last two terms in Equation (1) are of the same order of 
magnitude they tend to cancel out and thus Ah is also a simple function 
of the Froude number. This may explain the lack of variation in the 
backwater ratio with F, but it also raises a question. 

That question is —if a A 3 is large compared to the approach 


velocity head term in Eq. (1) (and it is not so in the model tests) will 
the backwater ratio still apply? 

Figures 11 and 12 can be readily combined for a given value of n into 
a single graph in which the backwater ratio is plotted against the con- 
traction ratio for equal values of the variable C,. The experimental 
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data, however, is quite meager for bridge types other than the base 
condition as is indicated by the scatter of data in Fig. 12. 

It is hoped that the authors will be perfectly frank in disclosing any 
serious defects in the form of presentation of data illustrated in Fig. 14. 
Such a graph gives a direct answer for backwater knowing the stage- 

‘discharge relation at the site, the waterway area under the bridge for 
any stage, and the conveyance of the bridge section in relation to the 
conveyance of the approach section. It does not give a measure of the 
depth y, which is the depth in the eddy area at the downstream side of 
the abutment and usually less than the normal depth of the unconstricted 
stream. But that depth, although very important for indirect discharge 
measurements, appears to have no practical significance on bridge de- 
sign. The model tests indicate that water surface elevations at other 
points along Section 3 will usually be slightly higher than elevation h,. 

Any one experienced in measuring velocities and depths of flow at 
bridges will know that the bed as a whole (not considering localized 
scour around piers) may lower during high stages of flow if the stream 
bed is alluvial and especially if it is non-cohesive. This increase in 
area apparently would reduce the required velocity but what would be 
the effect on backwater? This is one of the many problems which still 
need to be investigated, and may justify use of admittedly approximate 
solutions of the backwater depth at the present juncture. 
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FIG.13 - BACKWATER DEPTH RELATIVE TO NORMAL DEPTH, VERTICAL FACE 
CONSTRICTION WITH SQUARE-EDGED ABUTMENTS, b= 2.95 FEET 
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FIG. 14 - BACKWATER DEPTH RELATIVE TO NORMAL DEPTH, VERTICAL FACE 
CONSTRICTION WITH SQUARE-EDGED ABUTMENTS, b= VARIABLE 
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DISCUSSION OF MORNING-GLORY SHAFT SPILLWAYS: 
PROTOTYPE BEHAVIOR 
PROCEEDINGS-SEPARATE NO. 431 


BERNARD DONCLON,'—The Morning Glory Spillway has been used 
at Kingsley Dam to discharge water for operation a few months of 
every year since 1947. The Morning Glory Spillway has never been 
operated as a flood control structure since the reservoir has never 
reached the stage where it was necessary to discharge excessive water. 

The statement in the Proceedings-Separate No. 431 to the effect that 
no debris or ice has passed through the Spillway, as the flow has 
always flowed under the tractor gates, is correct. No water has been 
discharged through the Morning Glory Spillway during extremely cold 
weather as the film of water between the tractor gates and the guides 
freezes the tractor gates to the guides, and they cannot be moved. The 
Morning Glory Spillway is not operated during the winter when an ice 
cap has formed on the reservoir. 

The water in the Morning Glory Spillway discharge conduit stands 
at a depth of approximately 20 feet when the Spillway is not in operation. 
A diver made a complete inspection of the portion of the conduit under 
water in September, 1949, and found the surface of the concrete conduit 
to be in prefect condition with no signs of erosion or cavitation. The 
surface of the concrete above the water shows no wear. 

The Morning Glory Spillway has only operated under a very small 
percentage of capacity. The tractor gates are all working perfectly. 
The Spillway is operating with very little noise and no noticeable 
vibration. 


1. Resident Engr. — Sup’t., Central Nebraska Public Power and Irrig. 
Dist., Hastings, Nebr. 
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DISCUSSION OF MORNING-GLORY SHAFT SPILLWAYS: 
DETERMINATION OF PRESSURE-CONTROLLED PROFILES 
PROCEEDINGS-SEPARATE NO. 432 


MAXWELL W. WHITE,’ MURRAY B. McPHERSON,? J.M. ASCE.— 
The characteristics of shaft or morning-glory spillways in important 
structures are determined either by a model test of the actual structure 
involved or by an application of principles arrived at by previous studies 
on model shaft spillways or by more fundamental studies on sharp- 
crested circular weirs. The scaling up of results from these model 
tests to prototype conditions requires that the “scale effects” (i.e. . 
factors affecting flow over model weir which are not present in the 
prototype) either be negligible or else be allowed for in the calculation. 

The results given for discharge coefficients in Mr. Wagner’s paper 
cover spillways designed for maximum design heads from H,/D (where 


D = 2R) greater than 0.10 up past the submerged range. However, not 
many spillways have been designed to perform submerged. Of the 18 

morning-glory spillways for which details have been published® only 4 
are designed to perform submerged, 6 for ratios of H,/D between 0.10 


and 0.20 and 8 for ratios of H,/D less than 0.10. 


Tests were carried out on sharp-edged circular weirs in the Hydrau- 
lic Laboratory, Fritz Engineering Laboratory, Lehigh University be- 
tween 1949 and 1954.4 The project was sponsored by the Department of 
Civil Engineering and Mechanics. Professor W.J. Eney is Director of 
Fritz Laboratory, and Head of the Department. The tests were <imed 
at determining 


(1) the head-discharge relationship (i.e. between basic parameters 
Q, H and D); 

(2) the effect on this relationship of fluid properties such as surface 
tension and viscosity, particularly for small heads; 

(3) the point at which weir flow changes to flow through a reentrant 
tube. 


Research Assistant, Lehigh University, Bethlehem, Penna. 
Associate Professor of Civil Engineering, Lehigh University, 
Bethlehem, Penna. 

. “Morning-Glory Shaft Spillways: Prototype Behaviour,” by J. N. 
Bradley, Proceedings A.S.C.E. Separate No. 431, April, 1954. 
“Characteristics of Sharp-Crested Circular Weirs,” by M. W. 
White, G.M. Brey and R.G. Dittig. Report No. A-1342, Fritz 
Engineering Laboratory, Lehigh University, April 1954 (Available 
on loan). 


These results were then compared with those of other investigators 
to extend the range and value of the results. A brief survey of the 
investigations and of the results obtained, is reported here. 

The present results were not extended to pressure-controlled 
undernappes mentioned in Mr. Wagner’s paper, all tests being conducted 
with the undernappe aerated at atmospheric pressure and with the nappe 
springing clear from the crest. The approach conditions were also 
fixed at uniform radial flow. 

From preliminary tests, it was found, in agreement with the paper, 
that, as long as P was greater than D, the effects of the velocity of ap- 
proach head and the ratio H./ P were negligible on the discharge coef- 


ficient. Accordingly, this factor was eliminated from subsequent test- 
ing, in all cases P being greater than D. 

With these variables eliminated, the remaining factors can be 
grouped into 4 significant dimensionless parameters. 


Cr 


Np: a Reynolds No. 


N a Weber No. 


w’ 


go80 yi” 


Q 
The first parameter is a theoretical coefficient of discharge. 
Dimensional analysis shows that this coefficient is a function of H,/D, 


which expresses the contractive effect, a Reynolds No. which expresses 
the effects of viscosity and a Weber No. which expresses the effects of 
surface tension. 

In the tests at Lehigh University only one fluid (water) was used, 
therefore the last two parameters, Np and Ny» could not be varied 


independently, since both are a function of fluid properties and head, 
Hg, only. The effects of the two factors, then, could not be presented 
generally for any fluid nor could their individual effects be separated. 
For this reason the results presented here apply only to water, and, in 
the investigations, the two parameters, Ny and Np can be represented 
as a function of H, only, for water. 

Within the above limitations, the investigation fell into two parts, 
that is the effect on the coefficient of discharge, of 


(1) the head H, 
(2) the ratio H,/D 


i 

(2) Bs | 
D 

(3) go-s0 

(4) 

: 

f 

q 

i 

} 

| 
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The diameters tested at Lehigh University were 0.375 diameter, 
0.542' diameter, 0.708' diameter and 0. 875' diameter. Results were 
also obtained for weirs of 1. 13! diameter®; 1.14' diameter®; 1.66 
diameter’; 2.0' diameter, 3.97’ diameter and 5.96' diameter.® 

The original data of all calibration and experimental tests carried 
out at Lehigh University are on file in Fritz Laboratory. A summary of 
all data is available.‘ 

From all available experimental data, Fig. 17 was plotted. It shows 


the variation of apg with the head H., for various constant 


values of the ratio H,/D. For clarity, the ordinate has been plotted as 


Q ( Q 
), which, since curves are plotted 
D?-5g° 5 5p}. ¢° 5 D 


for constant values of Hs/D, has only the effect of separating the curves. 
The graph indicates a when H. is greater than approximately 0.10 
ft., the value of aD jazi constant for each value of H ;/D 


up to H./D ~0.30. That is, the'relationship is independent of H, (i.e. 


surface tension and/or viscosity) for heads greater than 0.10 ft., 
(corresponding to a Reynold’s No. of 16,600 and a Weber No. of 125). 

This provides more definite limits than Mr. Wagner’s paper on the 
minimum size of weir which may be used so that “scale effects” can be 
eliminated. The inconsistencies in the discharge coefficient for low 
heads, mentioned by Mr. Wagner on p. 16, were not found in the Lehigh 
tests, probably due to the smaller time required to establish steady 
conditions. The tendency of the discharge coefficient to “increase as 
head decreased” was due to surface tension and/or viscosity. The 
separate effects of surface tension and viscosity can be determined 
only by tests with different fluids flowing over the weir. Until such 
tests are;completed the present results apply only for water. 


Fig. 18 shows the relationship between arpa and H,/D for 
s 


the range where the former is independent of H,. On the log-log graph, 
this relationship plots as a straight line between the range of H,/D from 
0.03 to 0.30. The relationship is 


5. “Experiments on the Flow of Water Over Sharp-Edged Circular 
Weirs,” by H.J.F. Gourley, Proc. of the Institution of Civil 
Engineers, Vol. 184 Part I, 1910-11, pg. 297. 

. “Determination of Undernappe Over a Sharp Crested Weir, Circular 
in Plan, with Radial Approach,” by R. Dupont. Thesis submitted 
to Case School of Applied Science, Cleveland, Ohio, 1937. 

. Data received from Bureau of Reclamation, Denver, Colorado, of 
tests conducted by Mr. W.E. Wagner up till 29th Dec., 1951. 
“Tests on Circular Weirs,” by C. Camp and J.W. Howe, Civil 
Engineering, Vol. 9 No. 4, April 1939, pg. 247. (Original data kindly 
supplied by Professor J.W. Howe in January 1952—See Reference 4). 
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16005,0.80 


of Q 


This is in very close agreement with the values given by Mr. Wagner 
(pg. 13 and Fig. 9) for the more Limited range of his tests. 
The preceding equation can also be written 


Q = 1.66 Hi 


In Fig. 19, the parameter aq eps is plotted against the 
s . 


depth Hg. Once H, exceeds approximately 0.10 ft., the value is constant 


at 1.66. For heads less than 0.10 ft., the available results for all the 
weirs fall on a smooth curve for heads as low as 0.02 ft. In addition to 
indicating that the effects of surface tension and/or viscosity are cor- 
rectly represented by Reynolds and Weber Nos. involving the head only 
(as written), the graph also shows the value of the discharge coefficient 
for any head below 0.10 ft., when water is the fluid medium. 

The point at which weir flow is “flooded out,” the regime changing 
to flow through a reentrant tube, was found in the Lehigh experiments, 
to occur when H,/D was approximately equal to 0.30, representing the 


upper limit, above which weir flow cannot exist. However, if a vortex 
is formed, this causes the flow to diverge from weir flow at an earlier 
point. This is probably the reason why Mr. Wagner and earlier investi- 
gators found the limit for “free” discharge to be lower than 0.30. The 
effect of vortex formation is to reduce the discharge for a given head,° 
as shown in Fig. 20, which is a discharge-head curve for an arbitrary 
one ft. diameter circular weir. If the vortex is eventually destroyed 
with the establishment of tube flow the vortex causes a gradual transi- 
tion between the two types of flow. When a vortex was suppressed by 
means of simple radial baffles, then weir flow continued up to the upper 
limit at a vlue of H,/D ~0.30. Most morning-glory spillways have 


radial baffles to prevent the formation of vortexes.* No quantitative 
results have been obtained for the effect of the vortex in terms of the 
vortex strength. 

The conclusions from the Lehigh tests may be summarized as 
follows: 


For the flow of water over a sharp-edged circular weir, with the 
following limitations 


9. “How the Vortex Affects Orifice Discharge,” by C.J. Posey and 
H. Hsu, Engineering News Record, March 9, 1950, pg. 30. 
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1, radial approach conditions (no vortex), 
2. P greater than weir diameter, 
3. aerated flow (i.e. undernappe at atmospheric pressure), 


the relationship between the discharge, head over weir and weir 
diameter (up to point of “flood-out”) can be expressed as 


H, -0.04 
Q K {75 ve D 


or ve D!-% 


The value of K is a function of the head over the weir H,.- If H, exceeds 
0.10 ft., K is constant, equal to 1.66. If H, is less than 0.10 ft. the 


value of K may be obtained from Fig. 19 for heads down to 0.02 ft. 

The upper limit at which weir flow changes in character to flow 
through a reentrant tube is given by H,/D = 0.30; the formation of a 
vortex causes the flow to deviate from true weir flow at a lower value 
of this ratio. 
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PROC EEDINGS-SEPARATES 


The technical papers published in the past year are presented below. Technical-division sponsorship is indicated by an 
abbreviation at the end of each Separate Number, the symbols referring to: Air Transport (AT), City Planning (CP), Con- 
struction (CO), Engineering Mechanics (EM), Highway (HW), Hydraulics (HY), Irrigation and Drainage (IR), Power (PO), 
Sanitary Engineering (SA), Soil Mechanics and Foundations (SM), Structural (ST), Surveying and Mapping (SU), and Water- 
ways (WW) divisions. For titles and order coupons, refer to the appropriate issue of “Civil Engineering” or write for a 
cumulative price list. . 


VOLUME 79 (1953) 


SEPTEMBER: 260(AT), 261(EM), 262(SM), 263(ST), 264(WW), 265(ST), 266(ST), 267(SA), 268(CO), 269(CO), 270(CO), 271(SU), 
272(SA), 273(PO), 274(HY), 275(WW), 276(HW), 277(SU), 278(SU), 279(SA), 280(IR), 281(EM), 282(SU), 283(SA), 284(SU), 
285(CP), 286(EM), 287(EM), 288(SA), 289(CO). 


OCTOBER:» 290(all Divs), 291(ST)*, 292(EM)*, 293(ST)®, 294(PO)®, 295(HY)*, 296(EM)*, 297(HY)*, 298(ST)*, 299(EM)®, 
300(EM)*, 301(SA)*, 302(SA)*, 303(SA)*, 304(CO)*, 305(SU)*, 306(ST)*, 307(SA)*, 308(PO)*, 309(SA)®, 310(SA)*, 311(SM)*, 
312(SA)*, 313(ST)*, 314(SA)*, 315(SM)*, 316(AT), 317(AT), 318(WW), 319(IR), 320(HW). 


NOVEMBER: 321(ST), 322(ST), 323(SM), 324(SM), 325(SM), 326(SM), 327(SM), 328(SM), 329(HW), 330(EM)*, 331(EM)*, 
332(EM)*, 333(EM)°, 334(EM), 335(SA), 336(SA), 337(SA), 338(SA), 339(SA), 340(SA), 341(SA), 342(CO), 343(ST), 344(ST), 
345(ST), 346(IR), 347(IR), 348(CO), 349(ST), 350(HW), 351(HW), 352(SA), 353(SU), 354(HY), 355(PO), 356(CO), 357(HW), 
358(HY). 


DECEMBER: 359(AT), 360(SM), 361(HY), 362(HY), 363(SM), 364(HY), 365(HY), 366(HY), 367(SU)°, 368(WW)°, 369(IR), 
370(AT)©, 371(SM)©, 372(CO)©, 373(ST)°, 37-\(EM)°, 375(EM), 376(EM), 377(SA)©, 378(PO)°. 
‘ 


VOLUME 80 (1954) 


JANUARY: 379(SM)°, 380(HY), 381(HY), 382(HY), 383(HY), 384(HY)°, 385(SM), 386(SM), 387(EM), 388(SA), 389(SU)°, 390(HY), 
391 (IR), 392(SA), 393(SU), 394(AT), 395(SA)°, 396(EM)°, 397(ST)°. 


FEBRUARY: 398(IR)4, 399(SA)%, 400(CO)*, 401(SM)°, 402(AT)*, 403(AT)4, 404(IR) 4 405(P0)4, 406(AT)4, 407(sU)4, 408(sU)4, 
409(ww)4, 410(AT)4, 411(8A)4, 412(PO)4, 


MARCH: 414(ww)4, 415(su)4, 416(sm)%, 417(M)4, 418(AT)4, 419(SA)4, 420(SA)4, 421(AT)4, 422(SA)4, 423(CP)4, 424(AT)4, 
425(SM)4, 426(IR)4, 427(ww)4, 


APRIL: 428(HY)°, 429(EM)°, 430(ST), 431(HY), 432(HY), 433(HY), 434(ST). 
MAY: 435(SM), 436(CP)°, 437(HY)°, 438(HY), 439(HY), 440(ST), 441(ST), 442(SA), 443(SA). 


JUNE: 444(SM)®, 445(SM)®, 446(ST)®, 447(ST)®, 448(ST)®, 449(ST)®, 450(ST)®, 451(ST)®, 452(SA)®, 453(SA)®, 454(SA)®, 
455(SA)®, 456(SM)®. 


JULY: 457(AT), 458(AT), 459(AT)©, 460(IR), 461(IR), 462(IR), 463(IR)©, 464(PO), 465(PO)°. 


AUGUST: 466(HY), 467(HY), 468(ST), 469(ST), 470(ST), 471(SA), 472(SA), 473(SA), 474(SA), 475(SM), 476(SM), 477(SM), 
478(SM)°, 479(HY)°, 480(ST)°, 481(SA)°, 482(HY), 483(HY). 


SEPTEMBER: 484(ST), 485(ST), 486(ST), 487(CP)°, 488(ST)°, 489(HY), 490(HY), 491(HY)°, 492(SA), 493(SA), 494(SA), 
495(SA), 496(SA), 497(SA), 498(SA), 499(HW), 500(HW), 501(HW)°, 502(WW), 503(WW), 504(WW)°, 505(CO), 506(CO)¢, 
507(CP), 508(CP), 509(CP), 510(CP), 511(CP). 


a. Presented at the New York (N.Y.) Convention of the Society in October, 1953. 

b. Beginning with “Proceedings-Separate No. 290,” published in October, 1953, an automatic distribution of papers was in- 
augurated, as outlined in “Civil Engineering,” June, 1953, page 66. 

c. Discussion of several papers, grouped by Divisions. 

d. Presented at the Atlanta (Ga.) Convention of the Society in February, 1954. 

e. Presented at the Atlantic City (N.J.) Convention in June, 1954. 
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